function F = cal_integral(const, q_indicator, ExportProbQ)

    density_omega = const.density_omega;
    E_zmQ = const.E_zmQ;
    E_zvQ = const.E_zvQ;
    E_zxQ = const.E_zxQ;

%% Calculate integrals

    E_zm0 = density_omega'*(E_zmQ.*(1-ExportProbQ).*q_indicator);
    E_zm1 = density_omega'*(E_zmQ.*ExportProbQ.*q_indicator);
    E_zv0 = density_omega'*(E_zvQ.*(1-ExportProbQ).*q_indicator);
    E_zv1 = density_omega'*(E_zvQ.*ExportProbQ.*q_indicator);
    E_zx0 = density_omega'*(E_zxQ.*(1-ExportProbQ).*q_indicator);
    E_zx1 = density_omega'*(E_zxQ.*ExportProbQ.*q_indicator);    
   
    
%% Return

    integ.E_zm0 = E_zm0;
    integ.E_zm1 = E_zm1;
    integ.E_zv0 = E_zv0;
    integ.E_zv1 = E_zv1;
    integ.E_zx0 = E_zx0;
    integ.E_zx1 = E_zx1;
    F = integ;

end